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ABSTRACT 

Radio relics detected in the outskirts of galaxy clusters are thought to trace radio- 
emitting relativistic electrons accelerated at cosmological shocks. In this study, using the 
cosmological hydrodynamic simulation data for the large-scale structure formation and 
adopting a diffusive shock acceleration (DSA) model for the production of cosmic-ray 
(CR) electrons, we construct mock radio and X-ray maps of simulated galaxy clusters 
that are projected in the sky plane. Various properties of shocks and radio relics, 
including the shock Mach number, radio spectral index and luminosity are extracted 
from the synthetic maps and compared with observations. A substantial fraction of 
radio and X-ray shocks identified in these maps involve multiple shock surfaces along 
line of sights (LoSs), and the morphology of shock distributions in the maps depends 
on the projection direction. Among multiple shocks in a given LoS, radio observations 
tend to pick up stronger shocks with flatter radio spectra, while X-ray observations 
preferentially select weaker shocks with larger kinetic energy flux. As a result, the 
shock Mach numbers and locations derived from radio and X-ray observations could 
differ from each other in some cases. We also hnd that the distributions of the spectral 
index and radio power of the synthetic radio relics are somewhat inconsistent with those 
of observed real relics; a bit more radio relics have been observed closer to the cluster 
core and with steeper spectral indices. We suggest the inconsistency could be explained, 
if very weak shocks with Mg < 2 accelerate CR electrons more efficiently, compared to 
the DSA model adopted here. 

Subject headings: acceleration of particles — cosmic rays — galaxies: clusters: intra¬ 
cluster medium — methods: numerical — shock waves 
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Introduction 


Since the discovery of a bow shock in the periphery of the bullet cluster IE 0657-558 (jMarkevitch et al. 


2002 1. it has been well established that shock waves exist in and around galaxy clusters. Using 
cosmological hydrodynamic simulations for the large-scale structure (LSS) formation of the Uni¬ 
verse, the origin and nature of shock waves in the intraclus ter medium (ICM) as well as in the 


Pfrommer et al. 

2006; 

Kang et al. 

2007; 

Skillman et al. 

2008; 

Hoeft et al. 

2008: 

Vazza et al. 


200 . 2 : 


200911 . 


These studies demonstrated that abundant shocks are produced by supersonic flow motions during 
the process of hierarchical clustering of nonlinear structures, and that they could be classified into 
two categories, according to their locations relative to the host structures. External shocks are 
formed at the outermost surfaces surrounding clusters, filaments, and sheets of galaxies by the 
accretion of cool (T ~ 10^ — 10^ K), tenuous gas in voids onto those nonlinear structures. Since 
the accretion velocity around clusters can be as high as Uacc ~ a few xlO^ km s“^ and the sound 
speed of accreting gas is ~ 5 — 15 km s“^, external shocks are strong with Mach number as large 
as Mg ~ 10^. On the other hand, internal shocks are produced inside the nonlinear structures, 
where the gas has been already heated to high temperature by previous episodes of shock passage, 
so their Mach number is typically low with Mg < 10. While most of internal shocks have Mg < 3, 
those with 2 < Mg < 4 play the most important role in dissipating the shock kinetic energy into 


heat in the ICM (e.g., IRvu et al 


2003 


Kang et al.l 120071 ). 


Internal shocks could be further classified by their origins into a number of types, which in fact 
may not be mutually exclusive. Turbulent shocks are induced by tur bulent flow motion s in the ICM 
and they are expected to have very small Mach numbers, Mg < 2 (jPorter et al.ll2015l ). Turbulent 
motions in the ICM could be generated by several dif ferent processes: major or minor mergers . 


AGN jets, galactic winds, galaxy wakes, and etc (e.g., ISubramanian et al. 


2006 


2012 


Rvu et al.l l2008l . 


Vazza et al.l l2012bl : iBrunetti fc .lonesi l2014l ) . Bow shocks found in the outskirts of merging 


clusters are commonly re ferred as merger shocks (iRoettiger et al 


1999 


Markevitch fc Vikhlinin 


2007 : ISkillman et al.ll2013l h In most cases merger shocks are weak with Mg < 3 le.g.. lGabici fc Blasi 


20031 1. Infall shocks form by infall of the war m-hot intergalactic medium (WHIM; T ~ 10^ — 10^ 
K) into the hot ICM along adjacent filaments (jHong et al.ll2014l l. They could have relatively high 
Mach numbers, reaching up to Mg ~ 10, and thus the ensuing cosmic-ray (CR) acceleration can be 
more efficient than in other types of shocks. 

Including t he shock in the bu l let cluster, a number of shocks have been foun d in X-ray ob¬ 


Russell et al. 

2010: 

Akamatsu et al. 

2012; 


servations, shocks are detected as sharp discontinuities in the temperature and surface brightness 
distributions, and their physical properties including the sonic Mach number, M^, are estimated 
from the ‘deprojected’ temperature and density jumps (IMarkevitch fc Vikhlininll2007l ) . Most of the 
shocks detected in X-ray observations are weak with Mg ~ 1.5 — 3. 


In addition, shocks in the ICM have been detected by radio observations, especially as the so- 
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called radio relics (see, e.g.. iFeretti et alJl2012l : iBriiggen et alJl2012l . for reviews). Radio relics are 
usually found in the cluster outskirts around the virial radius, r^h- The radio emission is understood 
as synchrotron radiation from CR electrons with the Lorentz fa ctor of ^ 10^ — 10^ that are 


belie v ed to be accelerated at t he shocks associated with them fe.g. lEnsslin et al.lll998l : lBagchi et al 


20061 : Ivan Weeren et ahl |2010|) . High energy nonthermal particles can be produced via diffusive 


shock acceleration (DSA) at astrophysical shocks, such as interp lanetary and supernova remnant 
shocks as well as cluster shocks in collisionless tenuous plasma (iBelll Il978l : iBlandford fc Ostriker 
19781 : iDrurvI Il983l b Moreover, it has been show n that turbulenc e flow motions in the ICM can 
produce magnetic fields of up to ~ /rG level (e.g., IRvu et al.1120081) . 


Since the radio-emitting electrons with the Lorentz factor of 7 e ~ 10^ — 10® would not advect nor 
diffuse in the ICM more than ~ 100 kpc away from the shock surface before they lose th e energy du e 
to radiative cooling via synchrotron emission and inverse Compton (IC) scattering (e.g.. lKangll201ll ). 
it is commonly thought that their acceleration sites are likely to be close to where the synchrotron 
emission is seen. So the physical properties of shocks in radio relics are inferred from observed 
quantities such as the injection spectral index a,t the shock edge, a\n\ = + 3)/2(M^ ~ 


the spatial profile of surface brightness (jDrurvIIl983l : lEnsslin et al 


1998 


Kang et al.ll2012l i. Note 


that the spectral index for the synchrotron spectrum integrated over the downstream behind a 


steady planar shock is expected to be a ainj -|- 0.5 (e.g. lKandl2015bl : iKang &: Rvull2015l L In most 


cases, the Mach number of shocks in radio rel i cs was found to be in th e rang e of ^ 1.5 —4.5 ( e.g., 


Clarke fc Ensslin 


2006 


Bonafede et al. 


20091 : 


van Weeren et al 


in tne 

2010, 

range 

2012; 

01 IVIa ~ i.C 

Stroe et al. 

)-4.0 ( 

2014hi 


Pre-acceleration of thermal electrons to suprathermal energies and subsequent injection into 
the Fermi first-order process at shock have been one of the outstanding problems in the DSA 
theory. It is thought that the particle injection and DSA acceleration might be very inefficient at 
weak shocks {Mg < 3) because of small density compression across shocks (e.g., iMalkov fc Drury 
200l|). Especially, in the so-called “thermal leakage” injection model, proton injection is expected 
to be strongly suppressed at weak shocks; the relative difference between postshock, proton thermal 
and flow speeds is gre ater and so it is l ess likely for postshock protons to recross the shock front 
at weaker shocks (e.g.. lKang et al.ll2nn2l ). Although postshock electrons move faster than protons, 
they are tied to magnetic field fluctuations more tightly because of smaller rigidities (i.e. pth,e = 
{nie/mp So it was speculated that the DSA of electrons at weak shocks would not be 

efficient either, for instance, not enough to explain the observed radio flux of spectacular giant 
radio relics such as the Sausage relic. 

A pre-exiting population of electrons with the Lorentz factor of 7 e ~ 10 — 100 w as suggested 


as a possible solution to the low el ectron inject i on pr oblem at weak cluster shocks (jKang et al 


2013; iPinzke et al.ll2013l ). Moreover, iKang et al.1 (j2014l ) has suggested that a K-like distribution of 
suprathermal electrons may exist in high beta (/3 = Pg/Ps ~ 100) ICM plasmas just like in solar 
winds, and facilitate the electron injection at weak shocks. R e cently , usin g Particle-i n -Cell (PIC) 
simulations of weak shocks in high beta plasmas, ICuo et al.1 (120141 ) and I Park et al.1 (120151 ) have 
shown that some of the incoming electrons gain energy via shock drift acceleration (SDA) and are 
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reflected specularly toward the upstream region. Those reflected particles can be scattered back 
to the shock surface by plasma waves excited in the upstream region, and then undergo multiple 
cycles of SDA, resulting in a power-law type suprathermal population up to 7 e ~ 100. These studies 
suggest a possibility that “self pre-acceleration” of thermal electrons to suprathermal energies via 
kinetic plasma processes at the shock itself might provide seed electrons enough to explain the 
observed flux level of bright radio relics. 


The morphology of radio relics is, in some cases, observed to be elongated and arc-like with 
a sharp edge on one side. And some of radio relics are found as a pair in the opposite side of 
clusters. So they ar e often interpreted as p r oducts of binary merge rs fe.g.. lEnsslin et al.Ml998l : 


Roettiger et al.l 1 19991 : Ivan Weeren et al.l l2010l : Ide Gasoerin et al.l l2014l l . Several numerical studies 


have suggest ed that merger shocks with sufficient amount of shock kine tic energy flux could produce 
radio relics (jNuza et al.ll2ni2l : IVazza et al.ll2ni2al : ISkillman et al.ll2ni.‘ll ). In these studies, synthetic 
radio maps were constructed by identifying shocks in simulated clusters and modeling CR elec¬ 
tron injection/acceleration and magnetic field strength. However, there remain a few issues to be 
resolved, before the simple picture of merger shocks being the origin of radio relics is accepted. 


The first issue concerns the frequency of observed radio relics. Although structure formation 
simulations have demonstrated that shocks are produced frequently during mergers and they should 
last for the cluster dynamical time of tdyn ~ 1 Gyr, only about 10% of X-ray luminous clusters 
host some radio relics, putative merger shocks, an d the fraction of merging clusters with giant 
radio relics is even much lower (jFeretti et al.l 120121 ). Recently, it has been suggested that ICM 
shocks may light up as radio relics only when they encounter fossil relativistic electrons that are 
left over from either a previous episode o f shock/turbulence acceleration or a radio jet from AGN 
(jShimwell et al.ll2015l : iKang fc Rvull2015l l. In such a scenario, only a small faction of ICM shocks 
become radio-emitting structures for a fraction of the dynamical time (< O.ltdyn ~ 100 Myr). So 
the rareness of radio relics could be explained. Of course, this model can be justified only if the 
injection and acceleration of electrons at weak shocks in the ICM is very inefficient, so that the 
re-acceleration of pre-existing CR electrons must be required for the birth of radio relics. 


The second issue involves the discrepancies in the shock properties inferred from radio and X- 
ray observations of a few radio relics. Among several dozens of observed radio reli cs, only a fraction 
of shocks associated with them also have been detected in X-ray observations (jNuza et al.l l2012l : 


Bonafede et al.ll2012l . and references therein). In the case of the Toothbrush relic in IRXS J0603.3, 
for exmaple, the ra dio index was measured to be Oinj ~ 0.6 —0.7, indicating a ‘radio Mach number’, 
Afradio ~ 3.3 — 4.6 (jvan Weeren et al.ll2012ll . But the temperature and density discontinuities in X- 
ray observations suggest an ‘X-ray Mach number’, Mx < 2, and the positi on of shock identified in 
X-ray observations is shif ted from that i n radio observa t ions b y > 200 kpc (jAkamatsu fc Kawahara 


20131 : lOgrean et al.l l2013l ) . In addition, iTrasatti et al.l (j2015l ) suggested that for the radio relic in 
A2256, if the observed index ~ 0.85 measured between 63 and 1360 MHz is interpreted as the 
injection index, the radio Mach number can be estimated to be Mradio ~ 2.6, while the temperature 
jump measured in X-ray observations implies the X-ray Mach number, Mx ~ 1.7. In the case of the 
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Sausage relic in CIZA J2242.8, Ivan Weeren et alJ (|2010l i used the observed radio spectral index near 
the edge (shock surface), ainj ~ 0.6, to obtain a radio-inferred shock Mach number, Mj-adio ~ 4.6. 
But X-ray-inferred Mach number rep orted later by Suzaku and Chandra observations tu rned out to 
be lower with Mx-ray ~ 2.54 — 3.15 (jAkamatsu &: Kawaharall2013l : iQgrean et al.ll2014l l. Recently, 
however, IStroe et al.l (|2014al l estimated a steeper value, Oinj ~ 0.77, by performing a spatially- 
resolved spectral fitting, implying a weaker shock with Mradio ~ 2.9 in good agreement with X-ray 
observations. 


In iHong et al.l (|2014l . hereafter, Paper I), we studied properties of shock waves in the outskirts 
of simulated galaxy clusters using sets of LSS formation simulations. We find that in addition 
to merger shocks, infall shocks are produced in ICMs by continuous infall of density clumps (i.e. 
minor mergers) along filaments of galaxies to hot ICMs. Unlike weak bow shocks {Mg < 2) driven 
by major mergers, infall shocks do not show pairing structures and have higher Mach numbers 
{Mg ~ 3 — 10). In a few cases (e.g., Coma 1253-1-275 and NGC 1265), observed radio relics 


are thought to be associa ted with infall shocks (jBrown fc Rudnick 
Qgrean fc Briiggenl 120131 ). 


2011 


Pfrommer fc Jones 


2011 


In this paper, we consider a scenario in which suprathermal electrons are injected directly at 
the shocks without a help of fossil CR electrons, and accelerated to radio-emitting energies. Then 
the frequency of radio relics depends on the shock statistics, while the radio luminosity of radio 
relics is determine d by the shock kinetic energy flux and the assumed DSA efficiency model (e.g. 
Kang fc R vdl2013l ). In this scenario the rareness of radio relics, compared to the frequency of the 
ICM shocks, can be controlled by adjustin g the models f or the DSA acceleration efficiency and the 
magnetic held amplihcation. For instance, IVazza et al.l (j2012al ) showed, using structure formation 
simulations and generating mock radio maps of simulated cluster samples, that radio emission tends 
to increase toward the cluster periphery and peak around 0.2 — 0.5rvir, mainly because the kinetic 
energy dissipated at shocks peaks around 0.2Rvir and the Mach number of shocks tend to increase 
toward the cluster outskirts. Such hndings can explain why radio relics are rarely found in the 
cluster central regions. 

We study the properties of radio and X-ray shocks in the simulated cluster sample of Paper I, 
based on a DSA model in which suprathermal electrons are injected and accelerated preferentially 
at shocks with higher Mach numbers. We hrst calculate synchrotr on emission at shock s by modeling 
primary CR electron population based on recent DSA sir nulations (Kang &: Rvull2013 l and magnetic 
field distribution based on a turbulent dynamo model (jRvu et al.ll2008l l. We produce mock radio 
and X-ray maps by projecting the synchrotron and bremsstrahlung emission, respectively. We then 
identify “synthetic radio relics” in the radio map and extract their properties, such as radio and 
X-ray shock Mach numbers, Mj-adio and Mx, and their locations in the map. We examine the radio 
and X-ray properties of relic shocks, and attempt to understand discrepancies derived from radio 
and X-ray observations in some radio relics. 


In Section [2l we present numerical details such as the calculation of synchrotron emission 
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at shock, the construction of mock radio and X-ray maps, and the extraction of radio and X- 
ray shock Mach numbers. In Section [3l we discuss how the two-dimensional (2D) projection of 
three-dimensional (3D) shock distributions affects the radio and X-ray observations of ICM shocks. 
In Section 01 we describe the properties of shocks in 2D projection as well as the properties for 
synthetic radio relics. Summary follows in Section [5l 


2. Numerical Details 
2.1. Clusters and Shocks 


The construction of galaxy cluster sample and the identification of shocks are described in detail 
in Paper I, so here they are summarized only briefly. The LSS formation simulations adopted a 
standard ACDM model with the following parameters: baryon density fraction Dbm = 0.044, dark 
matter density fraction Ddm = 0.236, cosmological constant fraction Da = 0.72, Hubble parameter 
h = i7o/(100km s“^Mpc“^) = 0.7, rms density fluctuation us = 0.82, and the primordial spectral 
index n = 0.96. An updated version of a particle-mesh/Eulerian hydrodynamic code described in 


Rvn et al.l (jl993l ) was used for the simulations. Three sets of simulations were performed: (1) 16 
different realizations of LSS formation in a comoving cubic box of L = 100/i“^Mpc with 1024^ 
uniform grid zones, (2) 16 different realizations in a comoving cubic box of L = 200/i“^Mpc with 
1024^ zones, and (3) one realizations in a comoving cubic box of L = 100/i“^Mpc with 2048^ zones. 
The sets of (1) and (2) were performed in an adiabatic (non-radiative) way, while (3) includes a 
mild feedback from star formation and cooling/heating processes. In Paper I, we find that, in all 
three types of simulations, the statistics of the physical properties of clusters and shocks are similar. 


In the simulation data, clusters are identified as regions around local peaks in the spatial 
distribution of X-ray emissivity. Around each peak, the X-ray emission-weighted mean temperature, 
7x,ch is obtained for the spherical volume within r < r 2 oo- Here, r 2 oo ~ l-3rvir is the radius within 
which the gas overdensity is 200 times the mean gas density of the universe. Those with ksTx^d > 2 
keV for 100/i“^Mpc box simulations and ksTx^d ^ 4 keV for 200/i“^Mpc box simulations are 
selected as synthetic clusters, resulting in a total of 228 clusters from the three sets of simulations. 
Here, ks is the Boltzmann constant 

Shocks (actually, shock zones) are identified by a set of criteria given in Paper I, and their sonic 
Mach numbers, Mg, are estimated from the temperature jump condition, T 2 /T 1 = (5M^ — + 3) 

/(16M|), where T is the gas temperature. Hereafter, the subscripts “1” and “2” indicate the 
preshock and postshock quantities, respectively. Only shocks with Mg > 1.5 are considered. Note 
that a shock surface normally consists of many of these shock zones. The shock speed and shock ki¬ 
netic energy flux at shock zone are, then, calculated as vi = Ms( 7 i-th,i// 0 i)^'^^ and /km = (l/ 2 )piuf, 
respectively, where Pth,i is the preshock thermal pressure and 7 = 5/3 is the gas adiabatic index. 
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The energy flux of CR protons can be estimated as 

fcR,p = V{Ms) ■ /kin = v{Ms) ■ {l/2)pivl, (1) 


where is the CR accelera tion efficiency via the DSA process. For the efficiency, the values 

presented in lKang fc Rvd (120181 ) (also shown in Figure 2 of Paper I) are adopted. 


In Paper I, we hnd that the morphology of shock surfaces is quite complex due to the dynamic 
history of clusters, and in general a connected shock surface can consist of a number of grid zones 
of different types of shocks including merger shocks (see Figure 7 in Paper I). 


2.2. Modeling of Magnetic Field Strength 


To calculate radio synchrotron emission at cluster shocks, we need to model the strength of 
magnetic helds as well as the energy spectrum of CR electrons. The ICM is observed to be perme¬ 
ated with mag netic fields of a few uG in the cluster core region and a few x 0.1 pG in the cluster 
periphery (e.g. iBonafede et al.ll201ll : iFeretti et al.ll2012l i. Suggested ideas for the generation and 
amplification of magnetic fields in the ICM include processes during primordial phase transitions, 
Biermann battery mechanism, plasma processes at collisionless shocks, different types of turbulence 
dynamo, and the ejectiori of galactic winds and AGN jets (e.g. iDolag et al.l 120081 : IRvu et al.l 120081 . 


2012 


Widrow et al.ll2012l ). 


Here we adopt the turbulent dynamo model of iRvu et al.l (120081 ) , which assumes that turbulent 
flow motions are induced via stretching and compression of the vorticity generated behind curved 
surfaces of shocks in clusters, and the magnetic helds are amplihed by the turbulent hows. Then 
the strength of the resulting magnetic helds can be modeled in terms of the number of local eddy 
turn-over by the following htting formula: 






— 0(t/teddy) 


0.04 • exp[(t/feddy - 4)/0.36] for t/feddy < 4 
(0.36/41) • (t/teddy - 4) -b 0.04 for t/feddy > 4 


( 2 ) 


where cturb is the turbulent energy density and teddy = 1/|V x F| is the reciprocal of the vorticity 
calculated from the local how speed. The htting function represents the fraction of the turbulent 
energy transferred to the magnetic energy v ia turbulent dynamo, and it is derived from a magneto¬ 


hydrodynamic simulation (IRvu et al.ll2008l i. This model predicts that the magnetic held strength 
reaches to a few /rC in the cluster center and decreases to ~ 0.1 pG toward the cluster outskirts. 


Carilli &: Tavlor 

o 

o 

Govoni & Feretti 

o 

o 

Bonafede et al. 

2011) 










































2.3. Modeling of CR Electron Spectrum 


For the energy spectrum of CR e lectrons, we first assume that the CR acceleration at shock is 
described by the test-particle model (IDrurvlIlQSSl h since most of the shocks found in clusters are 
weak (see Introduction). Then, the momentum distribution of CR protons at the shock position is 
described by the power-law form for p > Pmin, 


fpiP) = fpOXP ® 


( 3 ) 


where p = pj{mpc) is the dimensionless proton momentum, nip is the proton mass, and q = 
3 (t /((7 — 1) is the spectral index. The shock compression ratio is determined by the sonic Mach 
number as u = (7 -|- l)M ^/[(7 — 1)M^ -|- 2]. The normalization factor fpQ at each shock zone is 
determined by setting the CR energy flux through the shock zone as 


/cR,p = V 2 -nipC^ + 1 - 1 ^ fp{p)dp^ 


( 4 ) 


where V 2 = vi/a \s the postshock flow speed and /cR,p is given in Equation ([T|) . We do not consider 
the re-acceleration of pre-existing CRs in this work. 

Here, Pmin = Pm] is the injection momentum above which particles can participate in the DSA 
process. According to recent PIC and hybrid simulations, both protons and electrons initially 
can gain energies via SDA while confined betwe en the shock front and the preshock region by 


scatt erings due to self-generated upstream waves (jCuo et al.ll2014l : iPark et al.ll2015l : ICanrioli et al 


20151 . see Introduction). Particles should have the rigidity [R = pcje) large enough to cross the 
shock front, i.e., p\a] 3pth,p ISOpth.ej in order to take part in the full first-order Fermi process. 
Here, pth.p = ^/^^^^^-pkB^ is the most probable momentum of thermal protons of postshock gas with 
temperature T 2 , while the most probable momentum of thermal electrons is pth.e = (?R-e/''R'p)^^^Pth,p- 
In fact. Pin] is expected to depend on the shock Mach number, obliquity angle, and shock speed. 
Here, for the sake of simplicity, we set = 0.01 in all shock zones regardless of the shock speed 
and Mach number and the ICM temperature. 

Since the DSA process operates on CR protons and electrons with the same rigidity in the 
same manner, the momentum distribution of primary CR electrons at shocks should follow that o f 


CR protons, except less efficient injection and radiative cooling (e.g.,lKang 


2011 


Park et al.ll2015l ). 


The injection rate of CR electrons is expected be much lower than that of CR protons, because 
postshock thermal electrons need to be pre-accelerated from the thermal momentum to pin] in order 
to get injected into the DSA process. Since the electron pre-acceleration is not yet fully constrained 
by plasma physics (despite of recent PIC simulations, see Introduction), it is often parameterized 
by the CR electron-to-proton ratio, K^/p. Different t ypes of observations have indicated a wide 


range of ilie/p 10 — 10 ^ (e.g., ISchlickeised 


2002 


Morlino et al.ll2009l b So we adopt the CR 


electron momentum distribution of fe{p) = Kejp ' fpiP) for P ^ V mm with ^ O-Ol- Adopting a 


higher value of K^jp, as suggested in some recent PIC simulations (|Guo et al.ll2014l : ICaorioli et al 
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201a), will increase the amplitude of fe{p) and therefore the synchrotron radiation flux uniformly 


for all shocks. 


At the same time of gaining energy via the DSA process at shock, CR electrons lose energy 
via synchrotron emission and IC scattering. As a result s, an equilib rium momentum exists where 
the momentum gain is balanced by the radiative losses (|Kangll2nill ): 

1 1/2 

, (5) 


Peq 


micvi 


mp\J Ae^q/21 


Bi 


r 2 I n2 

-“eff,! ^ '°efr ,2 


where all the quantities except Pgq are expressed in cgs units and the electron equilibrium mo¬ 
mentum is normalized as Pgq = peqjmpC. Here, takes account for the 

loss due to IC scattering of cosmic background photons (Rcbr = 3.24(1 -|- z)^/xG) as well as 
the synchrotron loss. Then, the CR electron spectrum at the shock location has the form of 
fe{p) = K^Ip ■ fp{p) ■ exp [-(p/Peq)^] ■ 


As the primary CR electrons are advected downstream in the postshock region, they also 
lose energy with the cooling time Rad = (0.54Gyr)(Rgg 2/5/rG)“^p“^. So the cutoff momentum 
due to radiative losses decreases with the postshock distance (d) away from the shock front as 
Pgut (d) oc Peq/d. As a result, the downstream, integrated electron spectrum steepens by one power 
of p as fe{p) oc for p > where the ‘break’ momentum, 

Pb, « 0.54 "Vmi) . (6) 


1 Gyr / \ 5 pG J 

is the lowest momentum up to which the postshock CR electrons have cooled down. In our sim¬ 
ulation sample, a shock is defined as discontinuity within a grid zone of thickness Al = 48.8— 
195.3/i“^kpc. The advection time for CR electrons to pass such thickness is Rdv = Aljv 2 ~ 
0.1 — 0.4 Gyr for V 2 ~ 500 km s“^, which is of the order of the dynamical time of typical clusters or 
the merger time. For the CR electrons emitting at GHz in pG magnetic field (p ~ 7), the cooling 
time becomes shorter than the advection time (Rad(il) < Rdv), so they have cooled down before 
existing the shock zone. Hence, we assume that the volume-averaged spectrum of CR electrons 
within each shock zone has the following form: 


fe{p) = ^ " 


for p < Pbr 

Kejp ■ fp0PhvP~^'‘^^^ ■ exp [-(p/Peq)^] P > Phr 


(7) 


Note that we made several assumptions and simplifications such as the DSA efficiency in 
Equation ([T]), the magnetic field strength, B, in Equation ([2]), the electron to proton ratio, Re/p; 
and the volume-averaged electron spectrum in Equation ([7|). 

Modifications to models for those will lead to rescaling of synchrotron radiation flux. So these 
may affect quantitatively some of the shock properties in the flux-limited, mock sample of radio 
relics, because a radio relic will consist of multiple shocks in the projected maps (see Section 2.4). 
But their impacts on the estimates of derived radio Mach numbers would be marginal, so the main 
results of this work should remain mostly unaffected (see Section 2.5). 
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2.4. Mock Radio/X-ray Maps 


We calculate the radio synchrotron emission at shock zones, using the magnetic field strength 
given in Equation ([2|) and the CR electron spectrum given in Equation ([7]). For a single electron 
with p, the synchrotron power at frequency zz is given by Py^e{p ,G) = (\/3e^i?sin0)/(meC^) x 
F {y jvc)^ where Q is the angle between the electron momentum and magnetic field directions, = 
3[(p (rup/me))^ + l](ei?sin0)/(47rmec) is the characteristic frequency, and F{x) = x 
{K^j^{x) is modified Bessel function). Then, the synchrotron emissivity (in cgs units of erg cm~^ 
s“^ Hz“^ str“^) in shock zone can be estimated as the sum of for all CR electrons of momentum 
spectrum fe{p) in the zone, assuming random angular distribution (i.e., (sin0) = vr/d), 


V37re^B2 ^ / 

\Mp) 


fe{p)d?P ■ 


( 8 ) 


Note that jy is in fact the volume-averaged emissivity, because fe{p) in Equation ([7]) is averaged 
over the volume of the grid zone taking account for radiative cooling in the postshock region. 


The calculation of the bolometric X-ray emissivity due to thermal bremsstrahlung in the ICM 
is straightforward and can be done with 

jx ~ (1.2 X 10“^®erg cm“^ s“^ str“^) > (9) 

\mpj 

where the He abundance is assumed to be 7.9 % by number (i.e., the mass fractions of X = 0.76 
and Y = 0.24 ) and the Gaunt factor gff « 1.2 for T > lO^K is used. Note that jx is calculated at 
all grid zones, while jy only at shock zones. 


Mock radio and X-ray maps of each synthetic cluster are constructed by projecting the syn¬ 
chrotron emissivity, jy, and the bolometric X-ray emissivity, jx, along a depth of H = 4r2oo- The 
choice of the depth should not affect our results, since both the synchrotron and X-ray emissions 
beyond r 2 oo from the cluster center are in general negligible. All line-of-sights (LoSs) are assumed 
to be parallel to each other, since the angular sizes of observed clusters and radio relics are usually 
much smaller than one radian. Assuming that the system is optically thin, the synchrotron/X-ray 
intensities, or surface brightnesses, are calculated by integrating the synchrotron/X-ray emissivities 
along each LoS, i.e., ly = /^os F dl and Ix = jx dl. 


Since the spatial distribution of shock zones in clusters is quite complicated, they can be 
trans lated into radio structures of different m orphologies, depending on the projection direction 
(e.g., IVazza et al.ll2ni2al : ISkillman et al.l 120131 ). So we choose 24 equally-spaced viewing angles for 
projection, resulting in 24 realizations of radio/X-ray maps for each of 228 synthetic clusters. As a 
result, a total of 5472 radio/X-ray maps are generated and used for the statistics below. 


Radio telescopes have finite resolutions, so in practice they measure the surface brightness 
convolved with telescope beams. If a beam has effective angular area 0^, then the measured 
synchrotron flux within the beam is approximately Sy^^^^ ~ Iy9‘^{l + z)~^, where t'obs = ^/(l + z) is 
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the redshifted observation frequency. In this study, we adopt z = 0 and a beam of = 15'' x 15" 
for synthetic radio observation. The beam size is chosen to comparable to those of future radio 
surveys such as those of the Square Kilometre Array (SKA). 

Hereafter, we set i^obs = 1-4 GHz as the representative radio frequency, and present ji .4 and 
(Si ,4 at the frequency as radio quantities for synthetic observation. For the calculation of radio 
spectral index, two frequencies of 142 MHz and 1.4 GHz are used (see the next subsection). 


2.5. “Derived” Mach Numbers 

With shocks of different Mach numbers and kinetic energy fluxes existing in clusters as men¬ 
tioned in Introduction, the radio/X-ray structures projected to the sky may consist of shock surfaces 
of different characteristics. So it may not be straightforward to define the properties for observed 
structures such as radio relics. In our 3D simulation data cube, each ‘shock zone’ is specified by 
the sonic Mach number, GR spectrum, and radio emissivity, while each zone is assigned with the 
gas density, velocity, temperature, magnetic field strength and X-ray emissivity, as described in 
previous subsections. In observations (real or synthetic), on the other hand, the physical properties 
of shocks must be extracted from 2D projections of radio/X-ray emissivities. Here, we describe 
the quantification of sonic Mach number Mg of shocks associated with 2D projected radio/X-ray 
structures. 

We specify the derived Mach numbers in two different ways. First, the weighted Mach numbers, 
M "4 and M^, are defined as the averages along LoS, weighted by their synchrotron and X-ray 
emissivities, i.e., = J^^gjiAMgdl/ /l^s ji. 4 d^ and = J^^^jxMgdl/fi^^gjxdl, respectively. 

We consider that these weighted Mach numbers represent the true properties of shocks associated 
with 2D projected structures. 

Second, the observed Mach numbers, and on the other hand, are designed to 

mimic the Mach numbers inferred from radio/X-ray observations. In 2D projected map, a ‘radio 
pixel’ is the pixel that has at least one shock zone along its LoS and so has > 0. For each 
radio pixel, is calculated from the integrated spectral index of synchrotron flux, i.e., a = 

—dlnSu/dlni/ = + 1)/(M°4®^ — 1), obtained between zzi = 142 MHz and 1^2 = 1-4 GHz, 

following the same way that radio observers usually interpret their observed radio spectra (see 
Table 1). Note that the value of M°^ is not sensitive to the choice of low frequency with our model 
of GR electron spectrum. However, this practice is justified only when the break momentum of 
the volume-integrated electron spectrum is <C 1 — 10, that is, when the break frequency in the 
integrated radio spectrum is < 100 MHz. 

For the calculation of the X-ray emission-weighted temperature is calculated along each 

LoS as Tx = f-^^gjxTdl /in mock maps of clusters. Pixels are tagged as ‘X-ray shocks’, 
if IA log Tx I > 0.11, as in the shock identification scheme for 3D volume (see Paper I). For X- 
ray shocks, is estimated from the temperature jump, i.e. Tx, 2 /Tx,i = (5M^*’®^ — 1)(M^*’®^ -|- 
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3)/(16 The larger of the values of estimated along the two primary (x and y) directions 

on projected maps is assigned as the observed Mach number, i.e., = max(M^'’®, M^y ). Again, 

only X-ray shocks with > 1.5 are considered. Note that a pixel identified as X-ray shock may 
not be a ‘radio pixel’ and vice versa. In fact, only a small fraction (~ 7%) of radio pixels are 
identified as X-ray shocks. Although the fraction may differ in real observations, this suggests 
that X-ray observations could miss a substantial fraction of radio shocks. Since it is caused by 
the smoothing in projection, this problem may be difficult to be overcome, even if the angular 
resolution of X-ray observation is improved. 


3. Projection Effects on Radio/X-ray Map 


In this section, we examine how the projection affects the radio/X-ray observations of shocks 
such as the morphologies and derived shock parameters. As a representative example, we consider 
a cluster from 100/i“^Mpc box simulation with 2048^ grid zones. This cluster is identical to 
the one appeared in the left panels of Figures 1 and also in Figure 7 of Paper I. It has the X- 
ray emission-weighted temperature kBTx,c\ ~ 2.7 keV and r2oo ~ 1.81/i“^Mpc. The shock that 
produces the largest amount of CR protons in this cluster is an infall shock with Mg ~ 5 and 
/CR ~ 1-5 X 10'^'^ erg s“^(h“^Mpc)“^ (see Paper I). 

Figured] shows projected maps of the cluster, viewed from four different angles, zoomed around 
the area of r < r2oo (dashed circle), where r is the distance from the cluster center in the projected 
plane. Here, the synchrotron flux, 5i,4, is shown as contours of black solid lines, which is superposed 
with the bolometric X-ray surface brightness, Ixj in gray-scale. The lowest contour level for the 
synchrotron flux is £'1.4,min = 10“^ mJy/beam, which could be detected by future radio observatories 
such as the SKA. The projection direction of each map is given in terms of polar and azimuthal 
angles, 0 and 4>, respectively. 


The four maps in Figure [T] exhibit morphologies of radio structures, distinct from each other. 
For instance, the radio map in the upper-left panel seems to show a paired structures on the 
opposite side of the cluster, which could be interpreted as radio relics due to a pair of merger 
shocks. Here, we tag the radio pixels associated with the left structure as “L-relic” (red color), 
while those associated with the right structure as “R-relic” (green color). In other three panels, if 
the shock zone with the largest synchrotron contribution along a given LoS belong to either L-relic 
or R-relic, we color the corresponding radio pixel as red or green, respectively. In the upper-right 
and lower-left panels, elongated radio structures with length >1 — 2 /i“^Mpc are composed of both 
red (L-relic) and green (R-relic) pixels. This illustrates that even a single connected radio structure 
in the sky may consist of a number of disconnected shock surfaces in the real 3D volume. So it 
could be misleading, if we try to extract the nature of underlying shocks only from the morphology 
or shock parameters inferred from radio observations. The radio map in the lower-right panel, on 
the other hand, has a number of small structures, and the structure located near the center co uld 
be interpreted as so-called “radio mini-halo” (see, e.g.. iFeretti et al.l 120121 : iBriiggen et al.ll2ni2l . for 
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reviews). Here, we do not intend to address the natures of observed paired structures or halo-like 
structures in detail, but we just point that the morphology of observed radio structures critically 
depends on the projection of underlying 3D structures. 


In Figure El we compare the spatial distributions of T, p, and Mg in a slice (left panels) 
with those of Tx, M°^, and in 2D projection (right panels). The slice passes through 

0.37 /i“^Mpc ~ 0 . 2 r 2 oo away from the cluster center, and contains the most radio-luminous shock 
zones for the L and R-relics, which are projected at (—0.4/i“^Mpc, 1.0/i“^Mpc) and (0.4/i“^Mpc, 
0.1h“^Mpc) in the 2D maps, respectively. They are identified as infall shocks with Mg ~ 5.7 
for L-relic and Mg ~ 3.9 for R-relic. The right panels are from the same map as that for 
the upper-left panel of Figure (H and the contours are drawn with the level of 5 'i. 4 ^min- For 
only pixels with S' 1.4 > S'i. 4 ,min are plotted, while for Tx and only pixels with 

7x > 7x,min = 10“^°erg s“^ cm“^ str“^ are plotted. Note that the bolometric X-ray surface bright¬ 
ness i n outskirts near r r^nn ranges typically 1 — 10 X 10“®erg s“^ cm“^ str“^ in observed clusters 
(e.g., Ettori fc Balestra 2009). 


Comparison of Mg in the slice map (left-bottom) and Mf^ in the 2D map (right-middle) indi¬ 
cates that the Mach number and location of shocks with Mg > 3 agree well in the two distributions, 
especially for the shocks associated with L-relic. This is because the radio spectrum in projected 
maps is governed by one or a few radio-luminous shock zones with high Mach numbers in a given 
LoS. So AL° 4 ® derived from radio observations could be a good proxy for the Mach number of 
radio-luminous shocks, implying that the properties of such shocks could be extracted reasonably 
well from radio observations. 


On the other hand, there are noticeable differences between the distributions of T /Mg in the 
slice maps (left-top/bottom) and those of Tx/M^^ in the 2 D maps (right-top/bottom). As seen in 
top panels, the distribution of Tx is smoother than that of T, partly because the contribution from 
the X-ray bright ICM dominates over the emissions from the WHIM in surrounding hlaments, and 
also because the projection (integration along LoSs) inevitably irons out any sharp features. The 
temperature jump across a shock in the Tx map looks reduced. We hnd that M^^ derived from 
X-ray observations tends to underestimate the actual shock Mach number (see the next chapter). 


We further examine in detail the region near L-relic, where the slice maps indicate a shock of 
Mg ~ 5 — 6 formed by the infall of the WHIM along a filament. There is ~ 100 /i“^kpc of positional 
shift between the locations of shocks in the Mf^ (right-middle) and M^^ (right-bottom) maps. 
This is because Mf^ picks up the infall shock, while due to the dominant contribution of X-ray 
from the hot and dense ICM along the LoS, M^^ picks up a foreground shock with smaller Mach 
number formed in the ICM (sse below). This may explain spatial offs ets between the shock surfaces 


inferred from radio and X-ray observations in some radio relics (e.g.. lAkamatsu fc Kawahara 
Qgrean et al. 2013 1. 


2013 


Figure [3] shows the distributions of radio and X-ray emissivities, ji.4 and jx, as functions of Mg 
for shocks within r < r 2 oo in our simulated clusters. Both distributions have rather large spreads. 
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but the following points are clear, Firstly, jx decreases with increasing M^, while ji .4 increases with 
Mg, peaks at Mg ~ 5 and decreases slightly for higher Mg. These behaviors can be understood as 
follows. Weaker shocks which tend to form in hot and dense regions near the cluster core produce 
more jx, while shocks formed in cluster outskirts with Mg ~ several accelerate CRs and produce 
ji .4 most efficiently. Secondly, on average, jx varies over two orders of magnitude, while ji .4 varies 
over 10 orders of magnitude. Although both radio/X-ray surface brightnesses on projected maps 
are dominated by contributions from a small number of bright shock zones along each LoS, the 
tendency is much stronger in the case of radio. 


4. Shock Properties Derived from Radio/X-ray Maps 
4.1. Shock Pixels in Projected Maps 

In this section, we first examine the surface brightnesses and derived Mach numbers of pixels 
with shocks, i.e S'!, 4 , M 1 . 4 , Ix and Mx, obtained in projected radio/X-ray maps. As noted in 
Section 2.5, radio pixels which contain shocks in radio maps may not be identified as shocks in X- 
ray maps, and vice versa. So when we study a correlation between any two quantities, for instance, 
M° 4 ® versus we use a subset of pixels in which both quantities are specified. 

Figure m displays the relations between 5 i ,4 vs Ix, M° 4 ® vs 51 , 4 , and vs Ix for pixels 

within r < r 2 oo from the cluster center. The bottom panel shows that bright X-ray shocks with 
large Ix tend to be weak with small the brightest X-ray shocks have < 2. The middle 

panel, on the other hand, shows that the synchrotron flux £' 1,4 of radio pixels increases with 
up to ~ 4 and then decreases for larger M/ 4 ®. These behaviors should be guessed from those of jiA 
and jx shown in Figured As a result, the correlation between £ 1,4 and Ix turns out to be rather 
poor with wide variations in the top panel. What is clear is that the shocks which are brightest 
in radio are not necessarily brightest in X-ray. As a matter of fact, brightest X-ray shocks have 
modest radio brightnesses, and vice verse. 

Figure [5] displays the relations among mock-observed Mach numbers, M/ 4 ® and M^'’®, and 
weighted Mach numbers, and M^, for pixels within r < r 2 oo- As noted in the previous chapter, 
the X-ray Mach numbers tend to be smaller than the radio Mach numbers (upper panels). Again, 
this is mainly because X-ray observations incline to pick up weaker shocks than radio observations 
along given LoSs. The correlation between Mi .4 and Mx is rather poor; the Pearson’s linear 
correlation coefficient of M/ 4 ® and M^^® is r(logM/ 4 ® ,logM^'’®) = 0.11 (upper-left panel) and that 
of M /'4 and M^ is r(log M/’ 4 , log M^) = 0.22 (upper-right panel). So it could be misleading if 
M^*’® is guessed from M/ 4 ®, and vice versa. It is expected that Mobs and M^ for the same band 
are correlated better. We hnd that M/ 4 ® is often smaller than as shown in the lower-left 

panel, but the correlation between them is quite good with r(log M/ 4 ®, log M/( 4 ) = 0.77. The 
correlation between and M^ in the lower-right panel, on the other hand, is still poor with 

r(logM^'’® ,logM^) = 0.24. This demonstrates that the smoothing in projection would make the 
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estimation of shock properties harder in X-ray observations than in radio observations. 


In Figure m we examine the radial distributions of observable properties of shock pixels. Both 
and tend to increase with radius, and shows a larger variance than Shocks 

with X-ray Mach number > a few would be rare within r < r 200 ) while shocks with radio 

Mach number up to several could be found. Radio bright shocks, for instance, those with 

<S'i .4 > 'S'i, 4 ,min = 10“^ mjy/beam, are found mostly in r > 0.2r2oo- This suggests that future radio 
observations could detect many more shock structures in cluster outskirts. On the other hand. 
X-ray bright shocks are preferentially located and so found close the center. The bolometric X-ray 
surface brightness of X-ray shocks follows lx oc (r/ r 2 oo)~^ with /3 ~ 3 — 5 for r/ r^nn > 0-3, w hich 
is consistent with the observed X-ray profiles in cluster outskirts fe.g.. lEttori &: Balestrall2009l i. 


4.2. Radio Relics in Projected Maps 

We build up a catalog of “synthetic radio relics” by finding connected structures in our radio 
maps, which meet the following conditions. (1) All pixels within structures have 5i.4 > S'i. 4 ,min- (2) 
A structure has at least five pixels. Note that the area of a single pixel is AA = {L/Ng)"^ {Ng = 1024 
or 2048). So hve pixels corresponds to, for instance, Amin ~ 0.05(/i“^Mpc)^ in simulations of 
L = 100/i“^Mpc box with 1024^ grid zones. Using these conditions, we have 9583 radio relics 
samples: 4626 from simulations of L = 100/i“^Mpc with 1024^ grid zones, 4850 from simulations 
of L = 200/i“^Mpc with 1024^ grid zones, and 107 from simulation of L = 100/i“^Mpc with 2048^ 
grid zones. 

We then assign the following quantities to synthetic radio relics. The average distance from 
the cluster center, (R), is dehned by the synchrotron-weighted average of the distance of the pixels 
that belong to a radio relic, (R) = X^pmei ^‘S'i.4/X^pixeiwhere r is the distance of pixels from 
the cluster center. The derived Mach numbers of radio relics are defined by the synchrotron/X-ray 
weighted averages of the corresponding derived Mach numbers of pixels, and M^. 

For example, = EpixdMu 5l.4/Epto|Sl.4 and = Eptol W U/Epix.l'x- 

The number of pixels with could be smaller than that with since is specified 

only at the pixels identified as X-ray shocks. In fact, about 40% of sample radio relics contain 
no pixel with so the corresponding X-ray Mach number is not assigned to them. The 

synchrotron power at 1.4 GHz and bolometric X-ray luminosity of a radio relic are calculated with 
Pi.4 = 47rAAX)pixei'S'i.4 and Lx = dvrAA X)pixei respectively. 

Figure [7] displays the relations among the derived Mach numbers for sample radio relics. The 
relations follow and so look similar to those for pixels in projected maps shown in Figure [5l but 
the correlations are tighter with r(log , log = 0.65, r(log 

= 0.64, r(logM°^®j.giig , log M "4 = 0.99, and r-(log 0.38. Especially, the 

correlation between and ^.gjjg is very good, indicating that the Mach number estimated 

from radio spectral index would be a fair representation of suitably averaged Mach number of shocks 
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reij . This is partly 


associated with radio relics. Most of our sample radio relics have > 2.5 . 

because with the model CR acceleration efficiency we employ (IKang &: Rvull2013l h the amount of 
CR electrons emitting synchrotron is very small for weak shocks with Mg <2.5 (see Figure 3 of 
Paper I). If the CR electron acceleration at weak shocks would be more efficient than in our model, 
radio relics with smaller could be more common. As noted in Introduction, Merger shocks 

are expected to have mostly Mg <3. In fact, we hnd that our synthetic radio relics normally 
involve projections of multiple shocks along LoSs, resulting in different morphologies for different 
viewing angles (see Figure [T]). In ~ 40% of sample radio relics, the brightest pixels with largest 5 i ,4 
include infall shocks along LoSs, suggesting that infall shocks may account for some radio relics 
with flat spectra. This seems reasonable in the sense that infall shocks could be the major sources 
of CRs in clusters, as discussed in Paper I. 

Figure El displays the radial distributions of Pi, 4 , and Lx- Upper panels 

show that most of our synthetic radio relics are found at (R) > 0 . 2 r 2 oo- Although there are 
wide variations, and tend to increase toward outskirts, and on average the radio 

Mach number is larger than the X-ray Mach number. The lower-left panel shows that the radio 
power at a given radial distance can vary widely, yet most powerful radio relics are found at 
0.3 < (R)/ r 2 oo ^ 0.7. This results mainly from the combined effects of the radial distribution of 
shock kinetic energy flux and Mg, the strong dependence of CR accelera tion efficiency on M g, and 
the geometrical increase of relic surface area toward outskirts (see also IVazza et al.l l20I2al l . The 
lower-right panel shows a strong radial decline of the bolometric X-ray luminosity of sample radio 
relics (not including the contributions from the background ICM) toward outskirts. This reflects 
mostly the radial decreases of gas density and temperature; the effect of the increase of the Mach 
number and relic surface area is not large here. 


4.3. Comparison with Observational Data 


We compare our synthetic radio relics with observed radio relics. Table I lists the properties 
of the observed radio relics we use; here, Zd is the redshift of associated clusters, /lls is the 
largest linear scale on th e sky , and a is the integrat ed radio spectral index. They are chosen from 


van Weeren et al.l ([20091) and iBonafede et al.l (|20I2l l with the following criteria; (I) Zd < 0.2 (our 


radio relic sample is constructed at z = 0), (2) ^lls P with Amin = 0.05(/i“^Mpc)^, and (3) 
both Pi .4 and a available. Note that a in Table I was calculated between 74 MHz and 1.4 GHz. 
Although some observations suggested that the radio sp ectra of radio relics cou l d have curvatures 


and a may increase at high frequencies > 10 GHz (e.g.. Ivan Weeren et al.ll2009l : IStroe et al. 


2013 


Trasatti et al.ll2015l ). here we ignore the effect. 


Figure [9| shows Pi ,4 vs (P) and a for synthetic (colors) and observed (filled squares) radio 
relics. The distribution of Pi ,4 of observed radio relics seems to fall within the range predicted by 
our synthetic observation. On the other hand, there are noticeable differences in the distributions 
of (P) and a between synthetic and observed radio relics. First, four observed radio relics are found 
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with (R) < 0.43 Mpc among twelve, while the fraction of synthetic radio relics with (R) < 0.43 
Mpc is very small, less than ~ 1%. Second, while the integrated spectral index of synthetic radio 
relics lies 1.0 < a < 1.38 (corresponding to Mg > 2.5), some of observed radio relics have steeper 
spectral indices. For example, the radio relics found at Abell 548b B and Abell S753 have a > 2 


{Mg < 1.7), and the radio relic at Abell 2345W has a k, 1.5 {Mg ^ 2.1). The discrepancies could 
be partly because weak shocks close to cluster core might be under-represented in our structure 
formation simulations due to the limited spatial resolution and omission of some non-gravitational 
processes. But it could be also possible that weak shocks accelerate CR electrons more efficiently 


than we assume here, as hinted in recent PIC simulations (see, e.g. iGuo et al.l l2014l : iPark et al 


2015|). 


On the other hand, the re-acceleration of “pre-existing” CR electrons, which is not consid¬ 
ered in this study, could play an important role in the p roduction of CR electrons at weak shocks 


mp 

with Mg < 2 ie.g. iKang et al.ll2012l : iPinzke et al.ll2013l i. Especially, if fossil relativistic electrons 


I 


are present in the form of isolated clouds, instead of being spread throughout the general ICM 
space, the relative rm^eness of radio reli cs can be understood, as discussed in the Introduction 
(jShimwell et al.ll2015l : iKang &: Rvull2015l i. In such a scenario, ‘lighting-up’ of radio relics is gov¬ 
erned mainly by the occasional encounters of ICM shocks with fossil electron clouds, instead of 
the Mach-number-dependent DSA efficiency and the shock kinetic energy ffux. As a result, the 
frequency and radio luminosity of radio relics may depend less sensitively on the shock Mach num¬ 
ber, and much weaker shocks could be turned on as radio relics with steeper spectral index. These 
features are more consistent with the observational data in Figure 9. Moreover, only a small frac¬ 
tion of IC M shocks become r adio sources for the period much shorter than the cluster dynamical 
time sale (jKang fc Rvul 120151 ). So the model naturally explains why radio relics are rare, while 


struc t ure formation simu lations predict shocks should be ubiquitous in the ICM (e.g. iRvu et al 


20031 : IVazza et al.l l2012al ) . The exploration of this scenario requires extensions of the models we 


employ here, so beyond the scope of this paper. 


5. Summary 


The existence of shocks in clusters of galaxies has been established through observations of 
radio relics due to synchrotron emission from shock-accelerated CR electr ons (e.g.. 


2012) as well as through X-ray observations of shock discontinuities (e.g.. iMarkevitch fc Vikhlinin 


Bonafede et al 


20071 ). However, it is not clear whether the properties of shocks inferred from radio and X-ray 
observations represent their true nature, because we can only measure the quantities projected 
onto the sky, integrated along LoSs. In fact, for a few relics such as the Tooth brush relic and the 
radio relic in A2256, radio and X-r ay observations have reported i nconsistencies in quantities such as 


shock Mach number and position (jAkamatsu &: Kawaharall2013l : iQgrean et al. 
2nL5l b 


2013 


Trasatti et al 
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In this paper, we explored a scenario in which electrons freshly injected and accelerated to high 
energies via DSA process at the structure formation shocks produce radio relics. We constructed 
synthetic maps of galaxy clusters in radio and X-ray by employing the cosmological hydrodynamic 
simulation data reported earlier in Paper I. First, the volume-averaged synchrotron emissivity at 
shock z ones was calculated by adopting the CR electron acceleration efficiency based on a DSA 


mode l (jKang fc Rvul 120181 ) and magnetic fields based on a turbulent dynamo model (jRvu et al 


20081 ) . The bolometric X-ray emissivity at grid zones was also calculated using the gas density and 


temperature from simulations. Then, mock maps of the synchrotron and X-ray surface brightnesses 
in 2D projection were produced by integrating volume emissivities along LoSs for simulated clusters. 
In the synthetic maps, the properties of shocks and radio relics, such as the shock Mach number 
and location, were examined in detail. 


The main findings can be summarized as follows. 

1) In most cases, radio and X-ray shocks in 2D maps are the outcomes of projection of multiple 
shock surfaces along LoSs, because shocks are abundant in the ICM with a mean separation of shock 
surfaces, ~ l/i“^Mpc. As a result, the morphology of shock distributions in 2D maps depends on 
the projection direction; for the same cluster, very different morphologies may turn out for different 
viewing angles. 

2) Synchrotron emissivity depends sensitively on the shock Mach number, especially with our 
DSA model for the CR electron acceleration, while bremsstrahlung emissivity depends on the gas 
density and temperature. Hence, radio observations tend to pick up shocks of Ms ~ a few to several 
along a LoS, while X-ray observations preferentially select weaker shocks {Mg < 2) with high density 
and temperature. Consequently, the properties of shocks in 2D projected maps could be different 
in radio and X-ray observations. The shock Mach number estimated with X-ray observations tends 
to be smaller than that from radio observations, if a radio relic consists of multiple shocks along 
LoSs. In addition, the location of shocks from X-ray observations could be shifted with respect to 
that from radio observations. 


3) For radio relics, the Mach number estimated from the radio spectral index seems to be a 
fair representation of suitably averaged Mach number of shock surfaces associated with them. On 
the other hand, the discontinuities in X-ray temperature tend to be smeared due to projection 
effects including possible multiple shocks and multiple ICM components, resulting in the X-ray 
Mach number lower than the real Mach number of the associated shock. 


4) When the properties of our synthetic radio relics are compared with those of observed ra¬ 
dio relics, there are clear differences in the statistics of radial location and spectral index; more 
radio relics have been observed closer to the cluster core and with steeper spectral indices than 
our synthetic observation predicts. This discrepancy may imply that weak shocks in high beta 
plasmas in fact ycele rate CR electrons more efficiently than we model h ere (see, e.g . 


2014; 

Park et al. 

20151. 

Pinzke et al. ( 

20131. the 


2pi5|). Alternatively, as shown in the previous studies of lKang et al.l (j2012l ) and 


Guo et al 


([2013), the re-acceleration of pre-existing electrons in the ICM may enhance the pro- 
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duction of CR electrons at weak shocks with Mg < 2. Finally, we could also conjecture that radio 
relics might be activated only when shocks encounter clouds, which contain fos sil electrons either ac¬ 
celerated earlier by shocks/turbulence or left over from old radio jets (see, e.g. IShimwell et al.ll2015l : 


Kang &: Rvull2015l i. This model may explain how very weak shocks can turn on the synchrotron 


emission and why radio relics are rare relative to putative shocks in galaxy clusters. 
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Table 1. Properties of Observed Radio Relics 


Source Name 


PlA 

(R) 

Ills 

a 

Reference 



[10^4 W Hz-1] 

[kpc] 

[kpc] 




Abell 2256 

0.058 

3.95 

300 

1100 

1.2 

Abell 1240 N 

0.159 

0.427 

700 

650 

1.2 

Abell 1240 S 

0.159 

0.730 

1100 

1250 

1.3 

Abell 2345 E 

0.177 

2.62 

890 

1500 

1.3 

Abell 2345 W 

0.177 

2.83 

1000 

1150 

1.5 

Abell 3667 

0.056 

17.4 

1950 

1920 

1.1 

Abell 548b B 

0.042 

0.250 

430 

370 

> 2.0 

Coma 1253+275 

0.023 

0.284 

1940 

850 

1.18 

Abell 2163 

0.203 

2.23 

1550 

450 

1.02 

Abell S753 

0.014 

0.205 

410 

350 

2.0 

Abell 115 

0.197 

16.7 

1510 

1960 

1.1 

Abell 610 

0.095 

0.444 

310 

330 

1.4 


Clarke Sz Ensslin 12006') 
Bonafede et al. 120091 
Bonafede et al. 120091 
Bonafede et al. 120091 
Bonafede et al. 120091 
Bbttgering et al. 119971 
Feretti et al. 120061 
Giovannini et al. 11991 1 
Feretti et al. 120041 
Subramanvan et al. 120031 
Govoni et al. (2001) 
Giovannini fc Feretti 120001 
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Fig. 1.— Projected maps of a simulated cluster with ksTx^d ~ 2.7 keV at z = 0, viewed from four 
different directions with 6 and (p listed. Contours show the synchrotron flux at 1.4 GHz, Si, 4 , in 
logarithmic scale, from 10~^ to ImJy/beam with interval of 10®'^, where a beam of 6^ = Ih" x Ih" 
is adopted. Grey-scale maps show the bolometric X-ray surface brightness, Ix, in logarithmic 
scale, from 10“® (white) to 10“^ ergs s“^ cm“^ str“^ (black). Dashed circles show r 2 oo- Hed and 
green colors differentiate the shock membership associated with the left and right “radio relics”, 
respectively, in the upper-left panel. The red crosses mark the X-ray center of the cluster. 
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Fig. 2.— Left panels: Distributions of gas temperature T (top), gas density p (middle), and 
shock Mach number Mg (bottom) in a slice which passes through ~ 0.2r2oo away from the center 
of the same cluster as in Figure [TJ Right panels: Projected maps of X-ray emission-weighted 
temperature Tx (top), Mach number estimated from the radio spectral index, M° 4 ® (middle), and 
Mach number of X-ray shocks detected in the Tx map, (bottom), with the same viewing 

angle as the upper-left panel of Figured} The yellow (top panels) and black (middle and bottom 
panels) contours correspond to 5 i. 4 ^min = 10“^ mJy/beam. For Tx and the pixels with 

Ix > .^x,min = 10 “^° ergs s“^ cm“^ str“^ are shown, while for the pixels with 5 i .4 > 5*1 4 , min 

are shown. Dashed circles show r 2 oo- 
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Fig. 3 .— Synchrotron emissivity at 1 . 4 GHz, ji.45 (red) and bolometric X-ray emissivity, jx, (blue) 
vs Mach number Mg for shock zones within r < r2oo from the center in simulated clusters. Con¬ 
nected lines and error bars are averages and standard deviations of log ji,4 and logjx within Mach 
number bins of [log Mg, log Mg A log Mg]. 
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Fig. 4.— 5 i .4 vs Ix, Mff vs 5 i. 4 , and vs Ix for shock pixels within r < r 2 oo from the cluster 
center in projected maps. Here, Si ,4 is the synchrotron flux at 1.4 GHz, Ix is the bolometric X-ray 
surface brightness, is the Mach number estimated from the radio spectral index, and 

is the Mach number of X-ray shocks detected in the Tx map. Colors code the relative frequency, 
ranging from 10“^'^ (blue) to 1 (red) in the top and middle panels and from 10“^'® (blue) to 1 (red) 
in the bottom panel. 
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Fig. 5.— Relations among observed Mach numbers, Mff and , and weighted Mach numbers, 
and (see the text for definitions) for shock pixels within r < r 2 oo from the cluster center 
in projected maps. Colors code the relative frequency, ranging from 10“^'^ (blue) to 1 (red) in the 
upper-left and lower-right panels, from 10“^ (blue) to 1 (red) in the upper-right panel, and from 
10“^ (blue) to 1 (red) in the lower-left panel. Pearson’s linear correlation coefficients between two 
derived Mach numbers at each panel are 0.11 (upper-left), 0.22 (upper-right), 0.77 (lower-left), and 
0.24 (lower-right). Black lines show the perfection correlation. 
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Fig. 6 .— Radial distributions of Mf 4® , , S' 1.4 , and Ix for shock pixels within r < r 2 oo from 

the cluster center in projected maps. Colors code the relative frequency, ranging from 10“^ (blue) 
to 1 (red) in left panels, from 10 “^'^ (blue) to 1 (red) in the upper-right panel, and from 10 “^'^ 
(blue) to 1 (red) in the lower-right panel. 
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Fig. 7.— Relations among observed Mach numbers, and and weighted Mach 

numbers, and Lhe text for definitions) for synthetic radio relics found within 

^ < ^200 from the cluster center in projected maps. Colors code the relative frequency, ranging 
from 10~^ (blue) to 1 (red) in the upper panels, from 10“^ (blue) to 1 (red) in the lower-left panel, 
and from 10“^'^ (blue) to 1 (red) in the lower-right panel. Pearson’s linear correlation coefficients 
between two derived Mach numbers at each panel are 0.65 (upper-left), 0.64 (upper-right), 0.99 
(lower-left), and 0.38 (lower-right). Black lines show the perfection correlation. 






- 31 



0.1 0.1 

<R>/r2oo 


Fig. 8.— Radial distributions of ./bfi 4 ®reiic 5 -^x^reiici synchrotron power Pi, 4 , and bolometric X-ray 
luminosity Lx for synthetic radio relics found within r < r2oo from the cluster center in projected 
maps. Colors code the relative frequency, ranging from from 10“^ (blue) to 1 (red) in the upper 
panels, from 10“^ (blue) to 1 (red) in the lower-left panel, and from 10“^'^ (blue) to 1 (red) in the 
lower-right panel. 
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Fig. 9.— Synchrotron power Pi ,4 vs synchrotron-weighted average distance (R) and integrated 
spectral index a for synthetic radio relics found within r < r 2 oo from the cluster center in projected 
maps. Filled squares are data points for the observed radio relics listed in Table 1. Colors code the 
relative frequency, ranging from (blue) to 1 (red) in the left panel and from 10 ~^ (blue) to 

1 (red) in the right panel. 




